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ABSTRACT  OF  THESIS 


Applications  for  Three-Dimensional  Computer 
Graphic  Cloud  Representations 
Produced  from  Satellite  Imagery 

This  research  focuses  on  applications  for  three-dimensional  cloud 
representations  »(hich  are  displayed  above  a  topographic  surface.  The 
computer  graphic  cloud-topography  displays  are  produced  using  actual 
GOES  satellite  data,  Defense  Mapping  Agency  topography  data,  and 
rawinsonde  upper  air  data.  Data  from  these  three  sources  is  used  to 
generate  three-dimensional  models  *i<iich  mathematically  define  cloud 
cover  and  topography  for  a  given  area.  By  applying  computet  graphics 
techniques,  these  three-dimensional  models  can  be  displayed  as  con¬ 
tinuous  tone  images  or  perspective  line  drawings.  Satellite  data  and 
topogra;ihy  data,  displayed  in  this  manner,  produces  a  realistic  depic¬ 
tion  of  cloud  cover  displayed  above  a  topographic  surface.  As  for 
applications,  these  displays  can  be  used  for  cloud  studies,  briefing 
aids,  and  flight  simulations,  v, 

Arthur  C.  Meade 

Department  of  Atmospheric  Science 
Coloracto  State  university 
Fort  Collins,  Colorad'j  80523 
Spring,  J965 
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CHAPTER  1 


Introduction 

In  recent  years,  computer  driven  image  processing  systems  have 
been  developed  »(hich  combine  computer  technology  with  graphics  tech¬ 
nology  to  produce  systems  capable  of  displaying  a  wide  variety  of 
digital  Imagery.  As  a  result^  many  graphics  software  packages  have 
been  designed  to  support  these  systems.  Some  of  the  more  advanced 
graphics  software  packages  are  capable  of  displaying  three-dimensional 
(3-D),  mathematically  defined  mo(tels  as  continuous  tone  images  or 
perspective  line  drawings.  In  the  work  covered  by  this  thesis,  com¬ 
puter  graphics  technology  is  applied  to  produce  integrated  >;loud- 
topography  displays  from  actual  GOES  satellite  data  and  topography 
data.  In  order  to  do  this,  clouds  and  topography  have  to  be  mathema¬ 
tically  defined  in  three  dimensions. 

In  this  research  process,  two-dimensional  (2-D)  satellite  images 
and  2-D^ topographic  images  are  used  to  produce  3-0,  mathematically 
defined,  cloud- topography  models.  High  resolution  (512  x  512), 
topographic  images  are  produced  n^ich  represent  the  earth’s  terrain 
over  a  4.27  x  4.27  degree  area.  These  topographic  images,  mapped  in  a 
Mercator  projection,  can  be  produced  for  all  areas  in  the  continental 
United  States.  GOES  infrared  and  visible  Images,  remapped  into  a 


Mercator  project  ion,  are  produced  »<hi'ch  represent  cloud  covrer  data 
over  the'  same  4.27  x  4.27  degree  area.  These  Mercator  projection  ima¬ 
ges  accurately  overlay  each  other.  By  displaying  these  images  in  suc¬ 
cession  using  an  image  processing  system,  cloud  cover  and  cloud 
features  can  easily  be  associated  with  effects  from  uiderlying 
topographic  features.  More  importantly,  by  having  these  images 
accurately  overlay  each  other,  cloud  cover  information  can  be  obtained 
with  respect  to  topography.  Information,  contained  in  the  above. 
Mercator  projection  images,  plus  .rawinsonde  upper  air  data  is  used  to 
mathematically  define  3-D  cloud-topography  models  for  sub-areas  within 
ttie  4.27  X  4.27  degree  areas.  These  3-D  cloud- topography  models 
represent  the  cloud  cover  and  topography  within  the  given  ajfa-area  and 
can  be  displayed  using  MOVIE. BYU  software. 

M3VIE.BYU  is  a  graphics  software  package  which  was  developed  at 
Grigham  Young  University.  It  is  used  in,  this,  research  to  manipulate, 
rotate,  illuminate,  and  display  these  mathematic^  models  using  the 
VAX  11/780  computer  and  the  COMTAL  image  processing  system.  MDVIE 
displays  these  3-D  models  on  the  COMTAL  as  either  continuous  tone  ima¬ 
ges  or  perspective  line  drawings. 

Three-dimensional  cloud-topography  models,  displayed  by  MDVIE,  . 
create  a  realistic  depiction  of  cloud  cover  above  a  topographic  sur-' 
face.  Because  of  the  display's  enhanced  information  content  and 
visual  appeal,  this  thesis  focuses  on  applications  for  satellite  data 
displayed  in  this  manner.  Applications  described  in  this  thesis 
include:  cloud  studies,  briefing  aids,  and  flight  simulations.  The  . 
following  five  chapters  discuss  the  data,  recources,  software,  and 


methodology  used  to  create  the  3-D  cloud-topography  model.  Chapter  7 
describes  the  applications  for  these  models.  Included  are  examples 
v#)ich  illustrate  how  the  cloud-topography  model  can  be  used.  Finally, 
Chapter  8  contains  a  summary  of  work ' covered  by  this  research. 
Suggestions  for  improvement  are  included  along  with  suggestions  for 
future  work  in  this  area. 


mAPTER  2 


Data 

Three  data  sources  are  used  to  create  the  3*0  cloud-topography 
models.  These  includes  Defense  Mapping  Agency  topography  data, 
rawinsonde  upper  air  data,  and  GOES  satellite  data.  These  data  sour¬ 
ces  along  with  data  acquisition  procedures  are  described  in  the 
following  sections. 

A.  Topography  Data 

The  Defense  Mapping  Agency  (DMA)  read  elevation  data  from 
li250,000  scale  topographic  maps  for  every  30  seco  ids  of  latitude  and 
longitude  in  the  continental  United  States,  parts  of  sbuthern  Canada, 
and  parts  of  northern  Mexico.  Individual  elevations  within  this  data 
set  were  read  in  feet,  to  the  nearest  increment  of  a).  These  values 
represent  elevations  at  specific  geographic^  points  and  are  not 
average  values. 

The  magnetic  tape  containing  this  data  was  obtained  from  the 
National  0<<ter  for  Atmosi^ric  Research  (nCAR).  'me  ta(»  consists  of 
1099  tape  blocks,  each  pf  which,  contains  all  me  elevation  data  for  a 
1x1  degree  square.  Data  is  included  for  a  region  from  2^N  to  Sl^'N 
and  from  l^w  to  60ow.  The  tape  blodcs  are  org»iized  starting  in  me 
soumwest  corner  of  me  region,  scan  eastward  to  the  eastern  boundary. 
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then  1  degree  north  from  ^e  start  and  eastward  again.  A  total 
of  28  of  these  west-eas*’  scans  are  made.  Tape  blocks  «^ich  only  con¬ 
tain  sea  level  data  are  omitted.  Figure  2-1  shows  regions  »*nere 
topographv  data  is  available. 

Each  tape  block  contains  a  latitude  and  longtiude  followed  by 
14641  consecutive  elevations.  The  latitude  and  longitude  specify  the 
geographical  location  oi  the  1x1  degree  squares.  The  elevations 
contained  in  a  tape  block  are  organized  starting  in  the  southwest 
corner  of  the  1x1  degree  square,  scan  eastward  to  the  eastern  boun¬ 
dary  of  the  1x1  degree  square,  then  move  30  seconds  north  from  the 
start  and  eastward  again.  A  total  of  121  west-east  scans  are  made. 
Boundary  values  in  each  1  x  1  degree  square  are  repeated  in  adjacent 
1x1  degree  squares  (Joseph,  1982). 

B.  Rawinsonde  Data 

The  National  Weather  Service  (NWS)  maintains  a  network  of 
rawinsonde  stations  which  send  rawinsondes  aloft  twice  daily  at  0000^ 
GMT  and  1200  GMT,  These  rawinsondes  contain  an  instrument  package  and 
a  radio  transmitter.  The  instrument  package  contains  Instruments 
which  measure  temperature,  pressure,  and  hunidity  throughout  the  tro¬ 
posphere  and  lower  stratosphere.  As  the  rawinsonde  ascends  through 
the  atmosphere,  the  radio  transmits  these  measurements,  back  to  the 
station  where  the  measurements  are  recorded,  encoded,  transmitted 
to  users,  Rawinsonde  data,  collected  in  this  manner,  is  divitted  into 
two  types*  mandatory  level  data  and  significant  level  data. 

Mandatory  level  reports  a>ntain  i^servations  from  ten  specific 
atmospheric  pressure  levels*  1000,  850,  7CK3,  5(X3,  400,  300,  250,  200, 


Figure  Z-1 

Available  topographic  data 
(adapted  from  Joseph,  1982) 
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150,  and  100  mb.  For  each  level,  temperature,  wind  direction  and 
speed,  dew  point  depression,  and  height  of  the  pressure  surface  are 
reported.  In  addition,  surface  information  and  tropopause  information 
are  included  in  the  reports. 

Significant  level  reports  contain  data  for  any  level  where 
important  changes  in  the  temperature,  humidity,  and  wind  profiles 
occur.  This  report  contains  two  sections,  one  for  reporting  tem¬ 
perature  and  humidity,  and  the  other  for  reporting  wind  profiles. 

The  CSU  Atmospheric  Science  Department  receives  this  data  from 
the  FAA604  weather  circuit.  The  data  is  received  via  a  three-meter 
antenna,  .processed  by  the  VAX  11/780  computer,  and  stored  on  disk.  In 
this  research,  only  mandatory  level  data  is  used. 

C.  Satellite  Data 

The  National  Oceanographic  and  Atmospheric  Administration 
(NDAA)  .-operates  two  geostationary  meteorological  satellites  (DOES) 
v*iich  crbit  the  earth  at  altitude  of  35000  kilometers  above  the  ^ 
earth’s  equator.  Since  the  satellites  orbital  period  is  synchronous 
with  the  earth’s  rotation  rate,  the  satellites  remain  stationary  rela¬ 
tive  to  fixed  locations  on  earth.  One  satellite  is  positioned  near 
TSOW  longitude  and  given  the  designation  of  GDES-East,  The  other  is 
positioned  near  1350w  longitude  and  designated  as  GOES-Mest  (Epstein 
^  al.,  1584).  Currently,  COES-East  is  not  operational  and  OOES-west 
is  positioned  near  llOOw  longitude. 

The  instrument  flown  aboard  the  GOES  satellite  is  the  Visible 
Infrared  Spin  Scan  Radiometer  (VISSR).  Tf«  VISSR  consists  of  a  mirror 
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that  is  stepped  to  provide  north-south  viewing  of  the  earth  while 
satellite  rotation  provides  west -east  scanning.  The  mirror  is  stepped 
one  notch  following  a  west-east  scan.  A  total  of  1821  of  these  west- 
east  scans  make  ip  a  full  frame  image  of  the  earth.  Normally,  a 
complete  full  frame  image  is  scanned  every  30  minutes,  24  hours  a  day, 
seven  days  a  week. 

Radiation  received  by  the  mirror  is  reflected  into  a  16-inch 
diameter  telescope.  A  fiber  optics  bundle  is  used  to  couple  the 
telescope  to  eight  visible  band  detectors  sensitive  to'  visible 
radiation  (.34  -  .70  ym).  Germanium  relay  lenses  are  used  to  pass 
infrared  radiation  through  optical  filters,  providing  a  10.5  -  12.6  ym 
band  pass,  to  the  infrared  detector. 

The  visible  detector  has  a  nominal  resolution  of  about  one  kilo¬ 
meter  while  the  infrared  detector  has  a  resolution  of  only  about  eight 
kilometers.  The  infrared  data  is  oversatnpled  in  the  west -east  scan. 

As  a  resuU,  the  sampling  resolution  is  four  kilometers;  in  the  west- 
east  direction.  Therefore,  at  the  sub-satellite  point,  visible  data  . 
has  a  one  kilometer  ground  resolution  and  infrared  data  has  a  8  x  a 
kilof'^ter  ground  resolution  (Hord,  1982).  Ground  resolution  decre  ases 
away  from  the  sub-satellite  point  due  to  the  curvature  of  the  eartlu 

Output  from  the  infrared  and  visibile  detectors  is  digitized  on 
the  satellite  and  transmitted  to  a  satellite  ground  station  at 
Mallops  Island,  Virginia.  There,  navigation  parameters  are  added. 

The  data  is  formatted,  calibrated,  and  retransmitted  back  to  the 
satellite  for  relay  to  other  users. 

As  the  data  is  relayed,  the  CSU  Atmospheric  Science  Depaftmert's 
Direct  Satellite  Readout  Ground  Station  (DSRGS)  collects,  processes,  . 
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and  stores  this  data.  The  DRSGS  basically  consists  of  the  following 
componentst  two  dish  antennas  to  receive  the  data,  two  POP  11/34  com¬ 
puters  to  process  the  data,  and  three  disks  to  store  the  data.  A  ten- 
meter  dish  antenna  is  positioned  to  receive  data  from  the  GOES-West 
satellite  and  a  five-meter  dish  antenna  is  positioned  to  receive  data 
from  the  GOES-East  satellite.  In  normal  operational  mode,  only  one  DEC 
POP  11/34  computer  is  needed  to  process  this  data  since  GOES-East  data 
is  transmitted  on  the  hour  and  half  hour,  Wiile  the  GOES-West  data  is 
transmitted  at  15  minutes  past  the  hour  and  45  minutes  past  the  hour.' 
The  other  DEC  POP  11/34  is  used  as  a  backup  to  insure  data  collection 
continues  if  one  of  the  components  within  the  system  fails  or  needs 
maintenance.  Data  received  by  DRSGS  can  be  stored  on  disk,  9-track 
magnetic  tape,  or  transferred  to  the  IRIS  system  where  the  data  can  be 
displayed  (Laybe,  1983).  A  block  diagram  of  the  ORSGSARIS  system  is 
shown  in  Figure  2-2. 


1 


IRIS 


Figure  2*2 

Block  diagram  for  DRSGS/IRIS  system 


CHAPTER  3 

Interactive  Research  imaging  System  (IRIS) 

The  Interactive  Research  Imaging  System  provides  the  data  manipu¬ 
lation  capability  and  the  image  processing  capability  to  create  and 
display  3-D  cloud-topography  models.  ,  The  Interactive  Research  Imaging 
System  consists  of  two  main  components,  the  VAX  11/780  computer  and 
the  COMTAL  Vision  One/20  image  processing  system.  Since  these  two 
resources  were  used  extensively  in  the  research  process,  they  warrant 
further  explanation.  In  addition,  a  description  is  included  on  how 
these  resources  are  used  to  display  the  3-0  cloud-topography  models. 

A.  VAX  11/780  Computer 

At  the  heart  of  the  Interactive  Research  Imaging  System 
(IRIS)  is  the  DEC  VAX  11/780  computer. ,  This  mid-size  con^uter 
features  a  32-bit  processor  whic^  possesses  sufficient  computing  speed 
for  many  meterplogical  and  image  processing  applications.  Perhaps  the 
VAX  11/780' s  most  attractive  feature  is  the  VAX/VMS  Operating  System. 
This  operating  system  allows  very  large  prograns  to  run  efficiently  by 
using  a  virtual  memory  system,  with  this  system,  the  VAX  11/780  uses 
disk  storage  to  extend  main  memory  onto  disk.  To  the  user,  it  appears 
the  systems  main  memory  is  unlimited  in  size  (Staff,  Digital  Equipment 
Corporation,  1982). 
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The  CSU  Atmospheric  Science  Department's  VAX  11/7&0  is  configured 
with  3  megabytes  of  main  memory  along  with  1.1  gigabytes  &f  disk 
storage  capacity.  Two  DEC  TU77  tape  drives  provide  additional  storage 
capacity  in  the  form  of  a  9-track  magnetic  tape  and  provide  a  medium 
to  enter  data  onto  the  system.  A  communication  link  from  the  VAX 
11/780  to  the  DRSGS’s  POP  11/60  provides  a  means  by  which  satellite 
data,  collected  by  DRSGS,'  can  be  transferred  to  IRIS.  Another  link 
from  the  VAX  to  the  COMTAL  Vision  One/20  image  processing .system  pro¬ 
vides  a  means  to  display  data  received  from  the  DRSGS'.  One  additional 
link  from  the  VAX  to  the  Colowrite  C-4300  (Optronics)  provides  a  means 
to  obtain  hard  copy,  high  resolution,  color  transparencies  of  imagery 
displayed  on  the  COMTAL  monitor.  A  block  diagram  of  the  DRSGS/lRIS 
system  is  shown  in  Figure  2-2. 

The  CSU  VAX  11/780  supports  two  high  level  languages,  VAX  FORTRAN 
and  VAX  PASCAL,  and  one  assembly  level  language,  VAX-11  Macro.  All 
software  used  to  create  and  display  the  3-D  cloud-topography  model  is 
supported  by  VAX  FORTRAN. 

B.  COMTAL  Vision  One/20 

Another  key  component  of  the  IRIS  system  is  the  COMTAL  Vision 
One/20  image  processing  system.  This  system  produces  high  resolution 
(512  X  512)  ii^gery  over  a  full  range  of  brightness  values  in  shades 
of  ^ey  or  pseudocolor.  Although  the  COMTAL  Vision  bne/20  is  a  stan¬ 
dalone  system  which  performs  a  multitude  of  image  processing  opera¬ 
tions  by  itself,  it  can  be  interfaced  with  the  VAX  11/7CO  as  host. 

This  feature  gives  the  COMTAlAAX  configuration  exceptional  image  pro¬ 
cessing  and  display  capability. 


Digital  imagery  (i.e.  imagery  which  is  obtained  by  partitioning 
the  area  of  one  inrege  into  a  finite  two-dimensional  array  of  small, 
uniformly  shaped,  mutually  exclusive  regions  called  pixels)  can  be 
created  and  displayed  on  the  COMTAL  monitor  using  an  image  plane  (512 
X  512  X  8  bits).  By  assigning  8  bits  to  each  pixel  in  an  image  plane, 
256  (2®)  grey  shades  bra  possible.  These  images  can  be  displayed  on 
the  COMTAL  monitor  as  monochrome  images  or  as  pseudocolor , images. 

Monochrome  images,  defined  by  eight  pixels,  can  show  256  vels 
of  greyness.  Each  pixel  making  up  the  image  can  be  black  (pixel  value 
=  0),  lignite  (pixel  value  =  255),  or  any  level  of  greyness  in  between. 

Pseudocolor  images  are  obtained  by  arbitrarily  adding  color  to 
the  monochrome  images.  These  pseudocolor  images  are  produced  by  the 
COMTAL ’s  pseudocolor  processor.  This  processor  accesses  a  red,  green, 
and  blue  pseudocolor  memory  which,  when  activated,  produces  signals  to 
the  red,  green,  and  blue  guns  in  the  color  monitor.  Each  of  the 
pseudocolor  memories  contain  256  8  bit  locations  which  correspond  to  a 


level  of  greyness  in  the  monochrome  image.  Each  location  in  a  pseudo¬ 
color  memory,  contains,  a  value  from  0-25>  specifying  the  intensity  of 
the  color  produced  by  the  gun.  Mhen  th>  memories  are  activated  by  the 
pseudocolor  processor,  the  values  of  piftels  making  up  the  monochrome 
image  are  indexed  into  the  red,  green. 

Then,'  the  red,  green,  and  blue  guns  are 
intensity  creating  a  pseudocolor  image 
red,  green,  and  blue  phosphors  on  the  screen. 

A  digital  image  can  also  be  represented  by  a  graphics  plane  (512 
X  512  X  1  bits).  Since  each  pixel  is  specified  by  only  one  bit,  two 


land  blue  pseudocolor  memories, 
each  fired  at  tfie  specified 
by  exciting  a  combination  of 
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(2^)  grey  shades  are  possible.  A  pixel  represented  in  this  way 
appears  either  black  (pixel  value  =  1)  or  »(hite  (pixel  value  v  oj. 
Digital  images  which  are  defined  in  this  manner  are  termed  "OOMTAL 
graphics"  or  "graphics".  These  graphics  often  contain  line  drawings, 
grid  lines,  or  labels.  Since  the  CDMTAL  stores  a  color  attribute  for 
the  graphic  as  a  whole,  an  entire  graphic  can  be  set  to  any  of  16 
colors.  Another  COMTAL  feature  allows  the  user  to  overlay  a  graphic 
on  a  monochrome  or  pseudocolor  image.  This  capability  gives  the  user 
the  ability  to  label  and  grid  monod^rome  and  pseudocolor  Images; 
without  destroying  any  of  the  images  original  information  content 
(Staff,  COMTAL  Corporation,  1982). 

C.  Display  ffethods 

The  VAX/COMTAL  systems  are  used  to  display  the  2-D  Mercator 
projection  imagery  and  the  3-D  cloud-topography  models^  This  section 
describes  how  the  two  systems  are  used  to  produce  these  displays. 

Mercator  projection  images  are  created  on  the  VAX  11/780  computer 
and  displayed  on  the  COMTAL  as  a  monochrome  images.  In  this  thesis, 
these  images  often  appear  with  one  or  more  COMTAL  graphics  overlaid 
onto  the  monochrome  image.  These  graphics  contain  either  labels, 
state  boundaries,  or  other  representations. 

Cloud- topography  models  are  created  on  the  VAX  11/7B0  computer 
and  displayed  on  the  COMTAL  as  pseudocolor  images*  COMTAL  graphics,  or 
a  combination  of  both.  A  3-D  model  displayed  as  a  pseudocolor  image 
contains  a  continuous  tone  representation  of  the  model.  A  3-0  model 
displayed  as  a  COMTAL  graphic  contains  a  contour  line  respresentatidn 
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of  tne  man.  mst  of  too  cloua-topogiaphy  flispjays  appearing  in  this 
thesis  «re  Oispiayea  on  the  COMTPL  ponltot  as  a  pseuCocolor  image 
overlaid  by  one  or  more  COmAL  graphlos.  men  displayed  as  a  pseudo¬ 
color  image,  :.opo9raphic  bases  are  colored  green  and  cloud  represen¬ 
tations  are  colored  mite.  To  increase  the,  display’s  perspective, 
contour  lines,  mich  are  «ltten  tc  a  COMTAl  graphic,  are  overlaid 
onto  the  pseudocolor  image,  the  CQMtPL  graphic  overlaying  the  clouds 
is  colored  mite,  mile  the  comSL  graphic  overlaying  the  base  is 
colored  green,  this  display  method  Increases  the  display’s  perspec¬ 
tive  and  visual  appeal,  plus  adds  a  grid  to  the  topographic  base  for 
quantitative  scaling., 


CHAPTER  4 

Software 

An  extensive  amount  of  software  was  used  in  the  process  of 
creating  and  displaying  the  3-D  cloud-topography  models.  Software 
falls  into  three  main  groups t  IRIS  Software,  hCVIE.BVU  Graphics 
Software,  and  User  Software  «^ich  was  written  by  the  author.  The 
major  programs  within  each  group  are  discussed  with  respect  to  purpose 
and  function. 

A.  IRIS  Software 

IRIS  applications  programs  are  designed  to  process,  manipu¬ 
late,  and  display  satellite  data  received  by  the  Direct  Readout 
Satellite  Ground  Station  (DRSGS).  These  programs  are  executed  using 
the  Transportable  Applications  Executive  (AE).  TAE  is  a  set  cf  exe¬ 
cutive  programs  «^ich  interact  with  the  user,  to  manage  the  execution 
of  applications  programs.  The  main  advantage  of  TAE  is  that  it  provi¬ 
des  a  menu  driven,  user  friendly  means  to  run  applica^-ions  programs. 
TAE  is  also  interfaced  with  an  IRIS  data  base.  This  data  base  is  used 
by  the  applications  programs  to  store  satellite  image  navigation 
parameters,  image  file  attributes,  and  other  information  describing 
satellite  data  sets.  In  this  research,  IRIS  applications  programs  are 
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usefli  1)  to  create  digital  images  from  raw  satellite  data,  2)  to 
display  the  digital  images  on  the  COMTAL,  3)  to  create  IRIS  data  base 
entrys  defining  image  attributes,  and  5)  to  plot  geopolitical  boun¬ 
daries  on  CDMTAL  graphics.  Applications  programs  performing  these 
functions  are  described  below. 

CSUGOES 

Raw  GOES  satellite  data,  which  has  been  received  by  the 
DRSGS,  is  read  from  either  disk  or  magnetic  tape  and  processed  by  the 
CSUGOES  program.  This  p^gram  produces  image  files  which  can  be 
displayed  on  the  COMTAL  monitor.  With  each  image  file,  automatic 
navigation  parameters  are  also  read  into  the  IRIS  data  base. 

WINDOW 

Data  processed  by  CSUGOES  can  be  displayed  on  the 
COMTAL  by  the  WINDOW  program.  WINDOW  creates  a  sub-image  by 
extr^ting  a  subset  of  data  from  an  exisiting  image.  The  program 
modifies  the  apparent  resolution  of  the  data  to  make  the  sub-image  fit 
exactly  into  a  512  x  512  image  plane.  If  the  original  image  was 
snaller  than  512  x  512  pixels,  WIMX3W  can  be  used  to  expand  the  image 
to  fit  into  the  512  x  512  image  plane.  If  the  original  image  is 
larger  than  512  x  512  pixels,  WINDOW  shrinks  the  image  to  fit  into  the 
512  X  512  image  plane.  The  resulting  sub-image  is  displayed  on  the 
COMTAL  monitor. 


Satellite  data  which  has  been  processed  through  CSUGOES 
can  be  remapped  into  various  projections  using  the  REMAP  program.  In 
this  research,  REMAP  is  used  to  rrmap  satellite  data  into  Mercator 
projection  images  covering  a  4.27  x  4.27  degree  area.  These  remapped 
Mercator  projection  images  can  be  displayed  on  the  CDMTAL  monitor  by 
using  WINDOW. 

OEFSCT 

OEFSCT  is  used  to  define  output  image  file  attributes 
for  the  REMAP  program.  Parameters  entered  into  the  DEFSCT  program 
define  the  location  of  the  4.27  x  4.27  degree  area  for  the  output 
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image  file.  Other  parameters  define  the  projection  of  the  output 
file.  These  parameters  are  stored  in  the  IRIS  data  base. 

GPPLOT 

Images  which  are  displayed  on  the  COMTAL.  can  be  gridded 
by  plotting  geopolitical  boundaries  on  them  using  the  GPPLXIT  program. 
These  boundaries  are  usually  written  to  a  COMTAL  graphic. 

B.  MDVIE.BYU  Software 

kJVIE.BYU  is  a  general  puniose  graphics  package  used  to 
create,  manipulate,  and  display  3>0  mathematical  models.  The’ software 
package  consists  of  a  total  of  seven  prograos,  of  vhich,  only  DISPLAY 
was  used  in  this  research.  DISPLAY  is  the  largest  and  most  uphisti*' 
cated  program  in  the  package  and,  as  the  name  implies,  is  used  to 
display  the  3-0  cloud-topography  models. 


DISPLAY 


The  particular  version  of  DISPLAY  used  in  this  research  was 
obtained  from  the  CSU  Electrical  Engineering  Department.  This  version 
was  originally  configured  to  interface  with  the  RAMTEK  9050  graphics 
device.  Since  no  RAMTEK  devices  were  available  at  the  CSU  Atmospheric 
Science  Department  and  no  COMTAL  interfaces  were  commercially 
available,  DISPLAY  was  modified  to  interface  with  tiie  COMTAL. 

Three-dimensional,  mathematically  defined  models  can  be  manipu¬ 
lated  and  displayed  using  the  DISPLAY  program.  Models  are  made  up  of 
a  series  of  three  and  four  sided  polygonal  elements  which  are  used  to 
define  model  surfaces.  Each  element  consists  of  three  or  four  nodes 
which  represent  intersection  points  where  the  sides  of  the  elements 
meet.  To  define  a  model,  the  X,  Y,  and  2  coordinates  of  each  node  in 
the  model  must  be  defined  and  stored  in  a  coordinate  array.  To 
complete  the  model  description,’  the  exact  sequence  in  which  nodes  are 
to  be  connected  to  form  polygons  must  be  specified  and  stored  in  a 
connectivity  array.  A  model  can  further  be  broken  down  into  parts. 
Parts  may  be  defined  by  grouping  contiguous  polygonal  elements. 

These  groups  are  specfied  and  stored  in  a  parts  array.  ITie  coordinate 
array,  the  connectivity  array,  and  the  parts  array  completely  describe 
the  model  mathematically  and  are  stored  in  a  geometry  file. 

Models  described  in  the  above  manner  can  be  viewed  on  the  COMTAl 
monitor  using  DISPLAY,  inhen  displayed,  models  appear  as  either 
continuous  tone  images  or  perspective  line  drawings.  A  total  of  39 
interactive  commands  are  available  to  the  user  to  specify  how  the 
model  is  to  be  displayed.  The  SCOPE  command, is  nsed  to  specify  the 
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output  device  and  to  specify  n^ether  the  model  is  to  be  displayed  as  a 
continuous  tone  image  or  perspective  line  drawing.  The  PARTS  command 
allows  the  user  to  withhold  model  parts, from  the  display.  Other  com¬ 
mands  are  available  to  manipulate  the  model.  The  ROTATION, 
TRANSLATION,  FIELD,,  and  DISTANCE  commands' allow  the  user  to  view  the 
model  from  virtually  any  position..  Still  other  commands  allow  the 
selection  of  color  and  shading  for  continuous  tone  imagery.  The  COLOR 
command  is  used  to  specify  colors  by  part  along  with  background  color. 
Colors  are  specified  as  mixtures  of  the  three  primary  colors,  red, 
green,  and  blue.  Model  illumination  can  be  specified  by  the  LIGHT 
command.  Up  to  four  light  sources  can  be  arbitrarily  defined  in 
space.  The  VIEW  command  is  used  to  display  the  model  on  the  COMTAL 
color  monitor  (Christianson  ^d  Stephenson,  1981). 

C.  User  Software 

A  total  of  eight  user  programs  are  used  in  this  research. 
Seven  of  them  are  used  to.  create  the  3-0  cloud-topography  models. 

The  other  program  is  used  to  illustrate  an  application  for  the 
cloud- topography  model.  User  software  is  ccmpatable  with  both  the 
KIVIE.BYU  Graphics  Software  and  the  COMTAL  image  processing  system. 
These  eight  programs  are  described  briefly  below. 
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DMA512 

l)MA512  is  a  user  interactive  program  used  to  create  2-D, 
Mercator  topographic  images.  These  images  represent  the  earth's 
terrain  and  cover  a  4.27  x  4.27  degree  area.  The  exact  location  of 
the  area  covered  by  the  topographic  image  is  specified  by  entering  a 
reference  latitude  and  longitude  which  defines  the  northwest  corner  of 
the  4.27  X  4.27  degree  area.  After  the  reference  latitude  and  longi¬ 
tude  is  specified,  DMA512  extracts  the  appropriate  elevation  data  from 
the  DMA  data  tape,  sc.ales  the  data,  and  formats  the  data  as  a  512  x  512 
COMTAL  image.  This  CDMTAL  image  is  written  to  an  output  file. 

Scaling  parameters  are  also  written  to  this  file.,  Information  con¬ 
tained  in  this  image  file  is  later  used  by  another  program  to  mathema¬ 
tically  define  3-D  topography  models. 

SilFTR 

After  the  REMAP  program  is  uSed  to  create  a  Mercator 
satellite  image,  the  resulting  image  may  not  exactly  overlay  the 
Mercator  topographic  image.  If  a  better  fit  is  desired,  SHIFTR  can  be 
used  to  shift  or  translate  a  satellite  image  left  or  right,  up  or 
down,  a  specified  nunber  of  pixels.  The  left-right  shift  is  specified 
,  by  a  nuirber,  X.  A  positive  value  for  X  shifts  the  image  to  the  right 
X  pixels  «blle  a  negative  value  shifts  the  image  left.  Similarly, 
the  up-down  ^ift  is  specified  by  a  number,  Y.  A  positive  value  for  Y 
shifts  trie  image  up  Y  pixels  while  a  negative  value  ^ifts  the  image 
left.  The  resulting  image  is  shifted  ^e  specified  nunber  of  pixels. 


DEFINE 


DEFINE  is  a  user  interactive  program  vnhich  is  used  in  the 
process  of  generating  a  3-D  cloud-topography  model.  DEFINE  is  used: 

1)  to  display  Mercator  images  on  the  COMTAL  monitor,  2)  to  specify 
input  and  output  file  names  used  in  the  process,  and  3)  to  specify 
parameters  «(hich  define  a  sub-area  of  a  Mercator  image  to  be  trans¬ 
formed  into  three  dimensions.  DEFINE  writes  these  file  names  and 
parameters  to  an  output  file.  By  having  this  data  in  a  file,  it  can 
be  accessed  by  other  programs  »^iph  are  used  in  the  process.  This 
alleviates  the  need  to  enter  file  names  and  parameters  more  than  once. 

iNTCuD 

INTCLO  is  a  user  interactive  program  designed  to  produce 
cloud  layer  files.  These  files  contain  cloud  top  height  information 
and  cloud  base  information  which  defines  cloud  cover  for  the  specified 
sub-area  of  a  Mercator  satellite  image.  INTQJ)  is  used:  1)  to  deter¬ 
mine  cloudy,  and  clear  pixels  within  a  sub-area  of  a  Mercator  satellite 
image,  2)  to  determine  cloud  top  heights  for  cloudy  pixels,  and  3)  to 
determine  cloud  base,  heights  for  cloudy  pixels.  INTCLO  produces 
three  cloud  layer  files,  each  of  which,  represents  cloud  cover  for  a 
layer  in  the  atmosphere.  Information  contained  in  these  files  is 
later  used  by  another  program  to  mathematically  define  3-0  cloud 
models. 

TOPGEN 

TDPGEN  is  a  user  interactive  program  designed  to  produce 
3-0  topography  models.  These  models  mathematically  define  the 
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topography  for  a  sub-area  of  a  Mercator  topographic  image.  TOPGEN 
extracts  the  appropriate  elevation  data  from  the  topographic  image, 
scales  the  elevation  data,  and  produces  a  geometry  file  »^ich  contains 
a  topography  model.  This  geometry  file  is  «nritten  in  a  MDVIE.BYU 
readable  format  and  can  be  displayed  on  the  COMTAL  using  the  DISPLAY 
program. 

CLDGEN 

CLOGEN  is  a  program  designed  to  produce  3-D  cloud  models. 
These  models  mathematically  define  cloud  cover  for  a  sub-area  of  a 
Mercator  satellite  image.  CLDGEN  reads  the  three  cloud  layer  files 
produced  by  INTCLD.  For  each  cloud  contained  in  these  files,  the 
program  produces  a  geometry  file  which  contains  a  3-D  cloud  model. 

Up  to  26  geometry  files  can  be  produced.  These  files  are  written  in  a 
K3VIE.BYU  readable  format  and  can  be  displayed  on  the  COMTAL  using  the 
DISPLAY  program. 

ICRGER 

.  »€RGER  is  designed  to  produce  a  3-D  cloud-topography  model. 
This  model  mathematically  defines  cloud  cover  aid  topography  for  a 
sub-area  of  the  Mercator  images.  ICRGER  produces  this  model  by  com¬ 
bining  the  3-D  topography  model  (TOPGEN)  and  the  3-D  cloud  models 
(CLDGEN).  MERGER  writes  the  mathematical  description  of  the  3-0 
cloud-topogr^hy  model  to  a  geometry  file.  This  file  is  written  in  a 
I0VIE.6YU  readable  format  and  can  be  displayed  on  the  COMTAL  using  the 
DISPLAY  program. 
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FLTRAK 

FLTRAK  is  used  to  plot  a  fli^t  track  on  the  3-D  cloud- 
topography  model.'  These  flight  tracks  mathematically  define,  in  3-D, 
an  aircraft’s  path  from  take-off  to  landing.  This  program  is 
discussed  further  in  Chapter  7. 


CHAPTER  5 


2-0  Image  Generation 

The  process  of  creating  a  3-D  cloud-topography  model  begins  by 
generating  a  set  of  three  Mercator  projection  images.  These  three 
images  contain  topography  data  from  the  DMA  data  set,  infrared  data 
from  GOES  infrared  satellite  imagery,  and  visible  data  from  the  GOES 
visible  satellite  imagery.  Data  contained  in  these  images  is  used  to 
define  terrain  and  cloud  cover  for  a  4.27  x  4.27  degree  area  bf  the 
earth.  It  is  the  main  concern  here  that  the  three  images  accurately 
overlay  each  other  since  it  is  from  Wiese  images  that  the  3-D  models 
will  be  generated. 

At  this  point,  it  must  be  emphasized  that  the  visible  satellite 
image  is  not  absolutely  necessary  to  generate  3*D  cloud-topography 
models.  The  main  function  of  the  visible  image  is  to  aid  cloud 
interpretation  ^d  image  navigation.  Therefore,  the  process  of 
generating  a  3-0  model  will  work  equally  well  using  only  the 
topographic  image  and  the  infrared  image.  However,  without  the 
visible  image,  Wie  accuracy  of  cloud  interpretation  and  image  naviga¬ 
tion  may  suffer.  The  remaining  part  of  this  chapter  explains  how  the 
Mercator  images  are  generated  and  navigated. 
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A.  Mercator  Topographic  Images 

The  generation  process  begins  by  the  selection  of  a  4.27  x 
4.27  degree  area  to  be  transformed  into  a  2-D  Mercator  topographic 
image.  Once  the  area  of  interest  is  selected,  the  DMA512  program  is 
used  to  create  the  topographic  image.  To  specify  the  area  of 
interest,  the  latitude  and  longitude  of  the  northwest  corner  of  this 
area  is  entered  into  the  DMA512  program.  Once  the  area  is  defined,  a 
512  X  512  array  of  elevations  is  extracted  from  the  DMA  data  tape. 
These  elevations  are  converted  to  meters  and  a  scale  factor  is  calcu¬ 
lated  by  dividing  the  highest  elevation  in  the  array  by  255.  The 
scale  factor  is  used  to  scale  each  elevation  in  the  array  from  0  to 
255.  Once  scaling  is  complete,  the  data  is  formatted  as  a  COMTAL 
image  and  is  written  to  an  output  file.  The  scaling  factor  and  the 
reference  latitude  and  longitude  —  also  saved  in  the  file  for  later 
use. 

A  Mercator  topographic  image  is  shown  in  Figure  5-1.  This  image 
.spresents  the  terrain  for  an  area  covering  southeast  Wyoming,  south- 
test  Nebraska  and  northern  Colorado.  A  COMTAL  graphic  containing 
state  boundaries  has  been  overlaid  onto  the  image.  To  show  the  reso¬ 
lution  of  the  image,  several  geographical  features  have  been  iabeled. 

Each  elevation  in  the  image  is  scaled  0  to  255  and  is  represented 
by  a  pixel  ^ey  shade.  Data  is  scaled  such  that  elevation  is  propor¬ 
tional  to  grey  shade.  Lighter  grey  shades  represent  higher  elevations 
and  darker  grey  shades  represent  lower  elevations.  Each  grey,  shade 
increment  in  Figure  5-1  represents  a  17.68  meter  change  in  elevation. 
Individual  elevations  in  the  image  can  be  calculated  by  multiplying 
the  pixel  grey  shade  value  by  the  scale  factor  117.68  meters). 


Figure  5-1 

Mercator  topographic  image 
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Each  pixel  in  the  image  represents  an  elevation  »(hich  is 
separated  by  30  seconds  of  latitude  in  the  north-south  direction  and 
30  seconds  of  longitude  iri  the  east-west  direction.  Since  the  lati¬ 
tude  and  longitude  of  the  pixel  in  the  northwest  corner  of  the  image 
is  known,  the  latitude  and  longitude  of  each  pixel  in  the  image  can  be 
calculated  by  referencing  the  pixel's  X  and  Y  screen  coordinates.  The 
pixel  in  the  northwest  corner  of  the  image  is  referenced  by  screen 
coordinates  X  =  0  and  Y  =  0.  X  screen  coordinates  increase  to  the 
east  (right)  and  Y  screen  coordinates  increase  to  the  south  (down). 
Since  each  pixel  is  separated  by  30  seconds  of  latitude  in  the:  north- 
south  direction,  the  latitude  of  a  pixel  can  be  calailated  by 
multiplying  its  Y  screen  coordinate  by  30  seconds  of  latitude  and 
subtracting  the  resulting  number  from  the  reference  latitude. 

Likewise,  the  longitude  of  a  pixel  can  be  calculated  by  multiplying 
its  X  screen  coordinate  by  30  seconds  of  longitude  and  siijtracting  the 
resulting  number  from  the  reference  longitude. 

Since  the  elevation  and  geographical  location  can  be  calculated 
for  each  pixel  in  the  image,  the  ima^  contains  all  the  information 
needed  to  define  a  topographic  surface  in  three  dimensions.  Data  con¬ 
tained  in  this  image  will  later  be  used  to  define  a  topogrophic  base 
for  the  3-D'  cloud- topography  model. 

B.  Mercator  Satellite  Images 

Now  that  the  topographic  image  has  been  generated  for  a  4.27 
X  4.27  degree  area,  GOES  satellite  imagery  can  be  collected  for  this 
region.  Figures  5^2  and  3-3  show  infrared  and  visible  satellite 


Figure  5-2 

Infrared  GOES  satellite  image  for  2132  GMT, 
31  October  1984 


Figure  5-3 

Visible  GOES  satellite  image  for  2132  GMT, 
31  Ocotber  1984 
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imagery  for  2132  GMT,  31  October  84.  This  data  was  collected  by  the 
DRSGS,  processed  using  CSUGOES,  and  displayed  using  WirCOW.  This  ima¬ 
gery  covers  a  portion  of  the  western  United  Sates  including  the  area 
covered  by  the  topographic  image  shown  in  Figure  5-1. 

The  next  step  in  the  process  is  to  remap  the  data  contained  in 
these  images  (Figures  5-2  and  5-3)  into  a  Mercator  projection  that 
will  accurately  overlay  the  topographic  image  (Figure  5-1).  This  is 
accomplished  using  the  REMAP  program.  Figures  5-4  and  5-5  show  the 
resulting  infrared  and  visible  images  after  the  remapping  process. 

These  images  were  displayed  using  wlhCpw  and  gridded  with  state  boun¬ 
daries  using  the  GPPLOT  program.  By  comparing  the  topographic  image 
(Figuare  5-1)  with  the  visible  image  (Figure  5-5),  it  can  be  seen 
that  the  images  closely  overlay  each  other. 

To  show  how  closely  they  actually  do  match,  elevations  greater 
than  3569  meters  have  been  blacked  out  in  Figure  5-6.  The  black  area 
has  been  made  into  a  COMTAL  graphic  which  can  be  used  as  a  template  to 
overlay  onto  the  visible  image.  Mhen  this  template  is  overlaid  onto 
the  image,  the  orographlcally  induced  clouds  covering  the  mountain 
ranges  should  coincide  closely  with  the  template.  Figure  5-7  shows  the 
visible  image  after  the  template  is  overlaid.  From  Figure  5-7,  it  can 
be  seen,  that  the  Images  overlay  each  other  to  within  several  pixels 
after  the  remapping  process. 

Many  timee,  the  remapping  process  is  not  as  accurate  as  the  above 
exanple  shows.  Figure  5-8  is  a  visible  image  which  has  been  remapped 
into  a  Mercator  projection.  In  this  case,  mountain  ranges  are  plainly 
visible.  By  overlaying  the  template  onto  this  image.,  it  is  apparent 


Figure  5-4 

Mercator  infrared  image  for  2132  GMT 
31  October  1984 


Figure  5-5 

Mercator  visible  image  for  2132  GMT 
31  October  1984 


Figure  5-6 

Mercator  topographic  image  with  template  over laid 
representing  elevations  above  3571  meters. 


Figure  5-7 

Mercator  visible  image  with  3571  meter  template 
overlaid. 
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Figure  5-8 

Mercator  visible  image  for  1930  GMT,  23  January  85 
with  3571  meter  template  overlaid. 


Figure  5-9 

Mercator  visible  image  for  1930  GMT,  23  January  85 
showing  3571  meter  template  overlaid  after  correcting 
navigation  error. 
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that  the  image  does  not  overlay  the  topographic  image.  To  correct  the 
descrepancy,  the  ^^IFTR  program  can  be  used  to  insure  a  more  accurate 
fit.  Figure  5-9  shows  the  same  visible  image  with  the  black  template 
after  the  visible  image  has  been  shifted  two  pixels  to  the  left  and 
nineteen  pixels  down.  This  image  accurately  overlays  the  topographic 
image. 

Since  the  satellite  images  accurately  overlay  the  topographic 
image,  cloud  cover  data  can  be  obtained  with  respect  to  topography. 
Data  contained  in  the  Mercator  infrared  image  will  later  be  used  to 
define  the  cloud  representations  that  will  be  displayed  above  the 
topographic  base. 


.  CHAPTER  6 

3-D  Model  Generation 

The  previous  chapter  explained  how  a  set  of  Mercator  projection 
Images  is  produced.  This  chapter  explains  how  these  images  are  used 
to  produce  3-0  cloud-topography  models. 

,A.  Si)-Area  Selection 

The  generation  process  begins  with  the  selection  of  a  sub- 
area  that  is  to  be  transformed  Into  three  dimensions.'  This  sub-area 
can  be  defined  using  any  of  the  three  Mercator  projection  images.  If 
the  user  is  interested  in  a  cloud  formation,  the  sub-area  may  be 
defined  on  a  satellite  image.  If  the  user  is  interested  in  the  cloud 
cover  over  a  particular  area,  the  sub-area  may  be  defined  on  the  . 
topographic  image.  The  process  for  selecting  a  sub-area  is  the  same. 

A  sub-area  can  be  selected  for  essentially  any  area  covered  by  a 
Mercator  projection  image.  The  sub-area  can  include  any  16  x  16,  32  x 
32,  64  X  64,  128  x  128,  256  x  256,  or  512  x  512  pixel  area  within  the 
Mercator  image i  however,  the  sub-area  must  remain  within  the  bounds  of 
the  image.  The  location  of  this  area  can  be  defined  using  the 
program,  DEFINE,  which  first  displays  the  desired  images  on  the 
COMTAL.  The  location  of  the  sub-area  is  defined  by  specifying  an  X 
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and  Y  screen  coordinate  for  the  northwest  corner  of  the  sub-area  and  a 
side  length.  The  X  and  Y  screen  coordinates  can  be  determined  using 
the  cursor  control  device.  By  placing  the  cursor  over  the  northwest 
corner  of  the  desired  sub-area,  the  cursor  control  device  will  return 
the  X  and  Y  screen  coordinates  of  the  cursor  position.  The  side 
length  for  the  sub-area  is  specified  by  the  number  of  pixels  along  one 
side  of  the  sub-area.  Choices  are  16,.  32,  128,  256,  snd  512. 

Once  the  X  and  Y  screen  coordinates  and  the  side  length  are  entered 
into  DEFINE,  a  box  is  drawn  around  the  selected  sub-area.  If  the 
displayed  sub-area  is  not  acceptable,  another  one  can  be  defined. 
Figure  6-1  shows  a  sub-area  displayed  on  the  Mercator  infrared  image 
(Figure  5-4}.  In  this  case,  the  sub-area  encompasses  the  entire 
image.  Parameters  defining  the  sub-area  aret  X  screen  coordinate  » 

0,  Y  screen  coordinate  *  0,  and  side  length  *  512. 

Since  KJVIE.BYU  software  is  only  configured  to  display  sub-areas 
to  and  including  64  x  64  pixels,  any  sub-areas  which  are  larger 
than  64  X  64  pixels  must  be  reduced  to  a  smaller  array  sire.  This  is 
done  by  specifying  a  reduction  factor  which  defines  the  amount  by 
«^lch  the  original  array  is  to  be  reduced.  For  example,  in  Figure 
6-1,  the  defined  sub-area  is  512  x  512  pixels.  A  reduction  factor  of 
16  averages  each  16  x  16  square  in  the  sub-area  to  a  single  smoothed 
pixel.  The  result  is  that  the  entire  512  x  512  pixel  sub-area  is 
reduced  to  a  32  x  32  array  of  pixels.  If  the  defined  sub-area  is  32  x 
32  pixels  or  64  x  64  pixels,  a  reduction  factor  can  optionally  be 
selected.  The  final  result  of  data  reduction  is  that  the  data 


Figure  6-1 

Mercator  Infrared  iniage  for  2132  GMT,  31  October  84 
with  white  box  showing  defined  sub- area. 


contained  in  every  sub-area  is  reduced  to  a  16  x  16,  32  x  32,  or 
64  X  64  array  of  pixels. 

The  reduction  factor  is  specified  at  the  time  the  sub-area  is 
defined  using  the  D^INE  program.  The  DEFINE  program  allows  the  user 
to  only  enter  a  reduction  factor  that  results  in  a  16  k  16,  3?  x  32, 
or  64  X  64  array  of  pixels.  After  the  reduction  factor  is  entered, 
DEFINE  stores  the  X  and  Y  screen  coordinates,  the  side  length,  and 
reduction  factor  in  a  file  where  the  data  can  be  accessed  by  other 
programs  used  to  create  the  3-D  cloud-topography  model. 

B.  Cloud/k}  Cloud  Threshold  Selection 

After  the  sub-area  is  defined,  cloud  cover  for  this  area  is 
defined  by  determining  which  pixels  are  cloudy.  Cloud/no  cloud  deter¬ 
minations  are  bdsed  on  a  simple  threshold  technique.  This  technique 
involves  using  either  a  threshold  radiance  count  from  the  infrared 
image  or  a  threshold  brightness  count  from  the  visible  image  to  make 
the  cloud/no  cloud  determination.  Pixels  with  radiance  counts 
(bri^tness  counts)  greater  than  or  equal  to  the  threshold  are  con¬ 
sidered  cloudy.  Pixels  with  radiance  counts  (brightness  counts)  less 
than  the  threshold  are  considered  cloud  free. 

The  cloud/no  cloud  threshold  is  specified  by  using  the  INTCLD 
program.  INTCLD  beT'ns  by  letting  the  user  display  the  desired  images 
on  the  COMTAL.  One  Jie  Images  are  displayed,  the  user  inspects  the 
imagery  and  interprets  what  is  cloud  and  what  is  not.  The  cursor 
control  device  can  be  used  as  an  aid  to  determine  the  cloud/no  cloud 
threshold.  By  placing  the  cursor  over  the  threshold  radiance  grey 
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shade  Cbri^tness  grey  shade)' the  cursor' control  device  will  return 
the  value  of  the  radiance  count  Cbrightness  count)  at  the  cursor  posi¬ 
tion.  This  value  is  entered  into  the  INTCLD  program. 

After  the  cloud/no  cloud  threshold  is  entered,  INTCLD  extracts 
the  visible  and  infrared  data  from  the  sub-area  defined  on  the 
Mercator  satellite  images.  If  a  reduction  factor  grea^r  than  one  was 
specified,  the  data  within  ^e  defined  sub-area  must  be  reduced  in 
size.  Recall  that  the  reduction  factor  specifies  the  amount  of  by 
vi^iich  the  original  sub-area  is  to  be  reduced.  Thus,  if  a  reduction 
factor  of  16  is  specified,  every  16  x  16  square  in  the  sub-area  is 
averaged  to  a  single  smoothed  pixel.  For  satellite  data,  data  reduc¬ 
tion  is  done  by  averaging  the  individual  pixels  within  a  square  using 
the  following  scheme*  If  fifty  -percent  or  more  of  the  individual 
pixels  in  a  square  are  equal  to  or  greater  than  the  cloud/no  cloud 
threshold,  the  resulting  smoothed  pixel  is  consictered  cloudy.  The 
radiance  count  Cbrightness  count)  for  this  cloudy  pixel  is  calculated 
by  averaging  the  radiance  counts  (brightness  counts)  of  those  indivi¬ 
dual  pixels  equal  to  or  greater  than  the  doud/no  cloud  threshold.  If 
less  than  fifty  percent  of  the  individual  pixels  within  a  square  is 
equal  to  or  greater  than  the  cloud/no  cloud  threshold,  the  resulting 
smoothed  pixel  is  considered  clear.  The- radiance  count  (brightness 
count)  for  this  clear  pixel  is  calculated  by  averaging  all  of  the 
pixels  within  a  square. 

After  all  the  data  in  the  sub-area  is  reduced  in  the  above 
manner,  the  result  is  a  16  x  16,  32  x  32,  or  64  x  64  array  of  smoothed 
pixels.  This  array  can  be  displayed  on  the  CQMTAt  as  a  reduced  image. 
Figure  6-2  shows  the  reduced  infrared  image  for  the  sub^-area  shown  in 


Figure  6-2 

Reduced  infrared  image  for  2132  GMT,  31  October  84 


Figure  6-3 

Reduced  infrared  image  for  2132  GMT,  31  October  84 
with  cloud/no  cloud  threshold  applied. 
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Figure  6-1.  Individual  pixels  within  the  reduced  32  x  32  array  have 
been  duplicated  by  a  factor  of  8  (i.e.  each  pixel  in  the  32  x  32 
array  is  duplicated  to  an  8  x  8  pixel  array,  thus  producing  a  236  x 
256  pixel  image).  This  is  done  only  to  enlarge  the  viewing  area. 

The  final  step'  in  the  cloud/no  cloud  determination  process  is  to 
apply  the  cloud/no  cloud  threshold  to  the  reduced  infrared  image. 
Figure  6-3  shows  the  same  reduced  infrared  image  in  Figure  6-2  after 
the  cloud/no  cloud  threshold  is  applied.  In  Figure  6-3,  all  smoothed 
pixels  in  the  reduced  infrared  image  that  are  below  the  cloud/no  cloud, 
threshold  have  been  set  equal  to  0  and  appear  black.  Radiance  counts 
associated  Mith  each  cloudy  pixel  in  the  reduced  infrared  image  are 
later  related  to  temperature  and  are  used  to  determine  cloud  top 
heights. 

C.  Cloud  Top  and  Base  Height  Calculations 

The  next  step  in  the  process  of  creating  a  3-D  cloud- 
topography  model  is  to  determine  cloud  top  heights  and  cloud  base 
heights  for  each  of  the  cloudy  pixels  within  the  reduced  infrared 
image  (Figure  6-3).  This  is  done  by  relating  pixel  radiance  counts  to 
temperature.  The  resulting  temperatures  are  then  related  to  height 
using  rawinsonde  sounding  data. 

The  radiance  counts  associated  with  a  GOES  infrared  image 
correspond  to  a  defined  set  of  equal valent  black  body  temperatures. 
This  information  is  summarized  in  a  calibration  table  which  contains  a 
range  Of  radiance  counts  and  their  corresponding  black  body 


temperatures.  Table,  6-1  shows  the  standard  infrared  calibration  table 
used  in  this  research.  The  accuracy  of  these  temperatures  is  limited 
to  two  to  four  degrees  Kelvin  (K),  generally,  nearer  two  degrees  K  in 
winter  and  four, degrees  K  in  summer  (Hord,  19B2).  The  INTCLD  program 
uses  data  contained  in  Table  6-1  to  calculate  cloud  top  temperatures 
for  each  cloudy  pixel  in  the  reduced  infrared  image. 

After  all  the  cloud  top  temperatures  are  determined,  the  next 
step  is  to  relate  these  temperatures  to  heists  using  an  atmospheric 
profile  of  temperature  versus  height!  Atmospheric  profiles  of  tem¬ 
perature  versus  height  are  obtained  from  rawinsonde  upper  air  data. 
Here,  a  sounding  is  selected  that  is  both  near  in  proximity  to  the 
area  of  interest  and  close  to  the  valid  time  for  the  satellite  image. 
Once  the  sounding  is  selected,  it  is  retrieved  from  a  disk  file  on  the 
VAX  and  the  mandatory  level  data  is  manually  decoded.  Decoded  data 
is  entered  into  a  file  which  is  read  by  the  INTCLD  program.  INTCLD 
processes  this  sounding  data  to  eliminate  surface  inversions  and  other 
inversions  that  may  occur  in  the  troposphere.  Missing  levels  are  also 
.  eliminated  in  the  process.  The  resulting  sounding  has  a  negative  lapse 
rate  throughout  the  0epth  of  the  troposphere,  this  insures  that  only 
one  cloud  top  height  can  be  assigned  for  a  given  cloud  top  , 
temperature. 

Figure  6-4  shows  a  temperature  versus  hei^t  profile.  This  pro¬ 
file  was  selected  for  determining  cloud  top  heights  for  the  sub-area 
shown  in  Figure  6-1.  This  profile  was  created  from  a  Denver, 

Colorado  rawinsonde  sounding  valid  for  1200  GMT,  31  (totober  84. 

Figure  6-5  shows  the  same  Denver  sounding  after  having  been  processed 
by  INTCLD.  .Notice  that  the  surface  inversion  has  been  eliminated. 
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Figure  6-5 

Temperature  versus  height  profile  for  Denver, 
Colorado,  1200  GMT,  31  October  84  after  inversion 
has  been  eliminated. 


After  the  sounding  is  processed,  INTCLD  creates  a  table  iKhich 
contains  an  entry  for  each  temperature  interval  in  the  sounding.  Each 
entry  in  the  table  containsi  the  temperature  (Tg)  at  the  bottom  of  the 
layer,  the  temperature  at  the  top  of  the  layer,  the  height  (Hb) 
at  the  bottom  of  the  layer,  and  the  lapse  rate  for  the  entire  layer. 
Table  6-2  shows  the  table  created  from  the  sounding  data  shown  in 
Figure  6-5.  INTCLD  uses  data  in  this  table  to  calculate  cloud  top 
heights.  For  each  temperature  associated  with  a  cloudy  pixel,  INTCLD 
finds  the  Interval  «^ere  the  temperature  falls.  Using  the  lapse  rate 
for  this  layer,  INTCLD  interpolates  the  height  difference  from  the 
bottom  of  the  layer.  This  height  difference  is  added  to  the  height 
(H0)  at  the  bottom  of  the  layer.  The  resulting  height  is  assumed  to 
be  the  cloud  top  heicfit  for  the  cloudy  pixel.  Cloud  top  heights  are 
calculated  for  all  cloudy  pixels  in  the  reduced  infrared  image. 

After  all  cloud  top  heights  are  calculated,  the  minimum  cloud  top 
height  is  calculated.  This  minimum  cloud  top  hei^t  is  assumed  to  be 
the  maximum  allowable  cloud  base  height  for  all  the  clouds  within  the 
reduced  infrared  image.  This  height  is  displayed  on  the  user  terminal 
and  the  user  determines  a  reasonable  base  height  for  the  clouds  in  the 
sub-area.  This  value  is  entered  into  the  program.  As  long  as  the 
value  is  below  the  maximum  allowable  base  height  it  is  accepted  by 
INTCLD. 

Finally,  after  the  cloud  base  height  is  entered,  INTCLD  creates 
three  cloud  layer  files.  These  files  contain  ddta  which  define  cloud 
cover  and  contain  a  cloud  top  height  and  cloud  base  height  for  each 


Table  6-2 

Interval  table  produced  from  temperature 
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cloudy  pixel  in  the  infrared  image.  Information  contained  in  these 
cloud  layer  filt'^  is  later  used  to  create  3-D  cloud  representations 
that  will  be  displayed  above  the  topographic  base. 

D.  3-D  Topography  Models 

At  this  point,  a  3-D  topography  model  is  created  using  data 
contained  in  the  Mercator  topographic  image.  This  model  mathemati¬ 
cally  defines  the  base  for  the  3-0  cloud-topography  model  and  is 
created  using  the  TOPGEN  program. 

The  process  of  generating  a  3-D  topography  model  begins  by 
reading  the  data  contained  within  a  sub-area  of  the  Mercator 
topographic  image.  To  find  the  location  of  this  sub-area,  TOPGEN 
reads  the  file  which  was  created  by  a^e  DEFINE  program.  This  file 
contains  a  reduction  factor  and  the  parameters  which  define  the  sub¬ 
area.  Next,  each  pixel  value  contained  within  the  defined  sub-area  is 
read  and  converted  to  an  elevation  using  the  scale  factor  that  was 
calculated  when  the  topographic  image  wsis  created.  If  a  reduction 
factor  greater  than  one  was  specified,  ^e  elevations  contained  within 
the  sub-area  are  averaged.  The  final  result  is  that  a  16  x  16,  32  x 
32,  or  64  X  64  array  of  elevations  is  created  representing  the 
topography  for  the  specified  sub-area. 

To  define  the  topography  model,  elevation  data  contained  in  the 
above  array  must  be  scaled  horizontally  and  vertically.  Horizontal 
scaling  involves  transforming  the  location  of  each  elevation  in  the 
array  to  a  model  X  and  f  coordinate,  vertical  scaling  involves  trans¬ 
forming  each  elevation  in  the  array  tc  «  model  2  coordinate.  Model  X 
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and  Y  coordinates  a»e  calculated  as  follows t  One  is  subtracted  from 
each  absolute  array  location  and  mutiplied  by  a  resolution  factor. 
Thus,  a  resolution  factor  of  4  nould  scale  a  32  x  32  array  of  eleva¬ 
tions  from  0  to  124  in  each  direction.  The  elevation  in  the 
southwest  corner  of  the  suh-area  would  have  an  X  and  Y  coordinate 
corresponding  to  (0,  0).  The  next  elevation  to  the  east  would  have  an 
X  and  Y  coordinate  corresponding  to  (4,  0)  while  the  first  elevation 
to  the  north  would  have  an  X  and  Y  coordinate  corresponding  to  (0,  4) 
etc.  After  nodel  X  and  Y  coordinates  are  ca^lculated,  the  actual 
distances  respresented  by  each  side  of  the  nodel  are  calculated. 
North-south  and  east-west  distances  are  averaged  and  this  average 
distance  is  divided  by  the  nutiter  of  units  along  one  side  of  ttie 
model.  The  result  is  a  model  scale  factor  which  specifies  the  number 
of  meters  represented  by  one  model  unit.  Model  Z  coordinates  are 
calculated  by  dividing  each  elevation  in  the  array  by  the  model  scale 
factor.  After  all  the  model  Z  coordinates  are  calculated,  the 
resulting  nodel  X,  Y,  and  Z  coordinates  arc  stored  in  a  coordinate 
array. 

Each  model  X,  Y,  arid  Z  coordinate  is  referred  to  as  a  node.  To 
complete  the  description  of  the  topography  nodel,  the  exact  sequence 
for  connecting  each  node  is  specified.  This  sequence  is  stored  in  a 
connectivity  array.  The  result  is  that  th*  topographic  surface  is 
defined  by  a  series  of  four  sided  polygons. 

Usually,  the  horizontal  distance  represented  by  a  nodel  is  much 
larger  than  the  vertical  distance.  As  a  result,  a  scale  nodel  would 
appear  very  flat.  To  bring  out  topographic  features,  a  vertical  scale 


factor  can  be  specified.  This  vertical  scale  factor  is  used  to 
exaggerate  the  model’s  vertical  scale.  To  aid  the  selection  of  a  ver¬ 
tical  scale  factor,  TDPGEN  calculates  a  model  height  to  width  ratio 
and  displays  this  ratio  on  the  user  terminal.  If  the  user  decides  to 
exaggerate  the  vertical  scale,  a  vertical  scale  factor  is  entered. 
Then,  a  new  model  height  to  width  ratio  is  calculated  and  displayed. 

If  the  new  ratio  is  not  satisfactory,  a  new  scale  factor  can  be 
entered.  After  the  desired  height  to  width  ratio  is  obtained,  each 
model  Z  coordinate  in  the  coordinate  array  is  multiplied  by  this  fac¬ 
tor. 

Finally,  a  model  base  is  added  to  give  the  surface  a  solid 
appearance  and  a  parts  array  is  specified.  The  resulting  coordinate 
array,  connectivity  array,  and  parts  array  is  written  to  a  geometry 
file.  This  file  contains  the  complete  mathematical  description  for 
the  3-D  topography  model.  The  file  is  written  in  a  MDVIE.BYU  readable 
format  and  can  be  displayed  using  the  DISPLAY  program.  Another  file 
is  also  written  which  contains  scaling  data  used  to  create  the  model. 
Data  in  this  file  will  be  used  later  to  scale  3-D  cloud  models. 

The  3-D  topography  model  can  be  viewed  on  the  COMTAL  monitor 
using  the  DISPLAY  program.  Figure  6-6  shows  a  topography  model  as  it 
appears  on  the  CDMTAL  monitor.  The  model  is  displayed  as  a  continuous 
tone  image  superimposed  with  a  gr^ic  which  contains  model  contour 
lines.  No  rotations  have  been  applied  and  the  resulting  view  of  the 
model  is  from  the  top  of  the  model.  This  model  represents  the 
topography  for  the  entire  area  shown  in  Figure  5-1.  The  vertical 
scale  has  been  exaggerated  by  a  factor  of  15}  however,  with  this  view. 


Figure  6-6 

Top  view  of  3-D  topography  model. 
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topographic  features  are  not  readily  noticable.  Normally,  the  3-D 
topography  model  is  not  displayed  and  is  done  here  for  illustrative 
purposes. 


E.  3-D  Cloud  Models 

The  previous  section  explained  how  the  topogr^hic  base  is 
created  for  the  3-D  cloud-topography  model.  This  section  explains  how 
the  cloud  representations  are  created. 

The  cloud  representations  are  created  using  the  CLDGEN  program. 
CLOGEN  begins  by  reading  the  cloud  layer  files  which  were  produced  by 
the  INTCLD  program.  After  a  file  is  read,  a  *'cloud"  is  extracted  from 
the  cloud  layer  file.  Here,  a  "cloud"  is  defined  as  a  group  of  cloudy 
pixels  which  touch  each  other  in  either  the  east-west  direction  or 
north-south  direction.  Information  defining  a  cloud  is  transformed 
into  X,  Y,  and  Z  coordinates  using  scaling  infoniiation  which  was  pro¬ 
duced  by  TOPGEN.  Further  data  manipulation  produces  a  coordinate 
array,  connectivity  array,  and  parts  array.  These  arrays,  «^ibh 
mathematically  define  a  cloud,  are  written  to  a  geometry  file.  This 
process  continues  until  all  clouds,  have  been  extracted  from  the  cloud 
layer  file.  After  a  cloud  layer  file  is  depleted  of  clouds,  the 
second  and  third  cloud  layer  files  are  processed  in  the  same  manner. 
Final  output,  for  this  process  is  a  series  of  up  to  26  geometry  files 
each  of.  »4iich  contains  a  single  3-0  cloud  model. 

The  3-D  cloud  models  *hich  were  produced  for  the  data  shown  in 
Figure  6-3  can  be  viewed  on  the  COMTAL  monitor  using  the  DISPLAY 
program.  Figure  6-7  shows  the  3-D  cloud  models  as  they  ^pear  on  the 


Figure  6-7 

Top  ew  of  3-D  cloud  models  for  2132  GMT, 
31  October  84 
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COMTAL.  These  models  have  been  displayed  as  a  pseudocolor  image 
overlaid  by  a  graphic  »(hidi  contains  model  contour  lines.  Again,  no 
rotations  have  been  applied  and  the  resulting  view  is  from  the  top  of 
the  model.  The  vertical  scale  is  exaggerated  by  a  factor  of  15.  This 
figure  depicts  the  cloud  cover  over  the  entire  region  shown  in  Figure 
6-1.  Notice  the  similarity  between  Figure  6-7  and  Figure  6-3. 

Again,  these  3-D  models  are  not  normally  displayed  and  are  shown  here 
for  illustrative  purposes. 


F.  3-D  Cloud-Topography  Models 

The  final  step  in  the  process  is  to  combine  the  3-D 
topography  model  with  the  3-D  cloud  models  to  produce  the  3-D  cloud- 
topography  model.  Geometry  files  containing  these  models  are  ojtttined 
using  the  MERGER  program.  Final  ouput  from  this  program  is  a  geometry 
file  which  contains  the  3-D  cloud-topography  model. 


The  3-D  cloud-topography  model,  which  was  produced  by  combining 
the  cloud  models  shown  in  Figure  6-7  and  tf.e  topography  model  shown  in 
Figure  6-6,  can  be  vie'«d  on  the  COMTAL  monitor  using  the  DISPLAY 
program.  Figure  6-8  siwws  the  3-D  cloud-topography  model  as  it 
appears  on  the  CDMTAL.  Again,  no  rotations  have  been  applied.  Figure 
6-0  graphically  illustrates  how  the  topography  model  (Figure  6-6)  and 
the  cloud  models  (Figure  6-7)  are  combined.  Figure  6-8  represents  the 
topography  and  cloud  cover  for  the  areas  covered  by  Figures  5-1  and 
6-1,  respectively. 

Figure  6-9  is  the  same  cloud-topography  model  as  shown  in  Figure 
6r8  after  rotations  have  been  applied.  The  resulting  view  is  from 


Figure  6-8 

Top  view  of  3-D  cloud-topography  model  for  2132 
GMT,  31  October  84. 


the  southeast  looking  northwest.  A  portion  of  the  Rocky  Mountains  is 
visible  on  the  upper  left  hand  side  of  the  model.  This  view  snows  the 
vertical  structure  of  both  the  terrain  and  the  clouds. 


WAPTER  7 


Applications 

The  aim  of  the  previous  four  chapters  was  to  explain  data, 
resources,  software,  and  methodology  used  to  create  the  3-D  cloud- 
topography  model.  This  chapter  Illustrates  how  the  3-D  model  can  be 
used  for  cloud  studies,  briefing  aids,  and  flight  simulations. 

A.  Cloud  Studies 

GOES  satellite  data,  displayed  as  digital  imagery,  is  repre¬ 
sented  by  a  regularly  spaced  array  of  pixels.  Each  pixel  in  the  array 
is  assigned  a  specific  grey  shade  based  on  the  numerical  value  of  the 
pixel.  The  resulting  picture  or  image  is  flat  scene  representation 
that  varies  to  grey  shade  from  pixel  to  pixel.  Satellite  data  which 
is  displayed  in  this  manner  is  commonly  used  by  meteorologists  to 
study  vertical  cloud  structure  and  cloud  development. 

To  study  vertical  cloud  structure,  infrared  images  are  often 
enhanced.  Image  enhancement  involves  altering  the  image  in  some  way 
to  facilitate  the  interpretation  of  its  information  content.  One 
method  of  image  enhancement  involves  Increasing  or  decreasing  contrast 
over  a  range  of  pixel  values.  Using  this  method,  enhancement  curves 
are  defined  which  assign  lighter  or  darker  grey  shades  to  specified 


pixel  values.  Another  method  of  image  enhancement  involves  the  use  of 
color  to  show  detail  in  the  imagery.  Using  this  method,  pseudocolor 
enhancement  curves  are  defined  which  assign  mixtures  of  the  three  pri¬ 
mary  colors  (red,  green,  and  blue]  to  specified  pixel  values.  Both 
methods  are  used  effectively  to  show  vertical  cloud  structure. 

A  time  series  of  satellite  images  is  often  used  to  study  cloud 
development.  Normally,  a  time  series  consists  of  two  or  more 
satellite  images  separated  by  a  time  step  ranging  from  minutes  to 
hours.  Cloud  development  can  tr  studied  by  comparing  these  images  in 
the  time  series  and  observing  changes  in  U)e  cloud  field. 
Alternatively,  cloud  development  can  be  studied  by  creating  a  movie 
loop  from  the  time  series  of  satrllite  images.  Here,  satellite  images 
are  displayed  in  rapid  succession  on  an  image  processing  system.  The 
result  is  a  time  lapse  display  of  clouds  and  weather  systems  which 
produces  the  illusion  of  continuous  cloud  growth. 

The  same  concepts  can  be  applied  using  the  3-0  cloud-topography 
model.  Displays  are  produced  from  the  cloud-topograf^y  model  using 
K)VIE.BfU  software.  These  displays  can  be  used  to  study  vertical 
cloud  structure  and  cloud  development  in  a  sindlar  fashion  to  those 
methods  discussed  for  2-0  satellite  imagery. 

Three-dimensional  cloud- topoOraphy  models  can  be  used  to  study 
vertical  cloud  structure.  To  bring  out  detail,  the  model's  vertical 
scale  is  exaggerated.  Vertical  model  scaling  is  similar  to  image 
enhancement  in  the  sense  that  it  facilitates  the  interpretation  of  the 
model's  information  content.  By  exaggerating  the  vertical  scales 
small  height  differences  are  amplified  allowing  relatively  minor  cloud 


features  to  be  seen.  When  the  model  is  displayed,  these  cloud 
features  can  be  viewed  and  studied  from  virtually  any  angle  or 
perspective.  In  addition,  the  satellite  data  is  scaled  vertically 
proportional  to  actual  hei^t  rather  than  to  temperature. 

Three-dimensional  cloud-topography  models  can  be  used  to  study 
cloud  development.  Here,  a  time  series  of  cloud  models  is  used 
instead  of  a  time  series  of  satellite  images.  The  time  series  of 
models  can  be  displayed  on  the  CDMTAL  monitor  using  the  DISPLAY 
program.  Cloud  development  and  growth  can  be  studied  by  comparing  the 
models.  Alternatively,  cloud  development  can  be  observed  by  creating 
and  displaying  movie  loops  using  the  time  series  of  3-D  models. 

Figures  7-1,  7-3,  and  7-5  show  a  time  series  of  GOL'S  infrared 
images  as  they  appear  on  the  CDMTAL  monitor.  Individual  images  are 
separated  by  a  time  interval  of  approximately  15  minutest  The  series 
shows  the  development  of  cumulus  convective  activity  in  north-central 
Kansas  from  2215  GMT  to  2265  GMT  on  21  June  1986.  Figures  7-2,  7-6, 
and  7-6  show  the  corresponding  3-D  cloud-topography  displays  for  tha 
boxed  sub-area  on  the  2-D  infrared  images.  The  vertical  scale  has 
been  exaggerated  by  a  factor  of  five  to  show  the  vertical  cloud  struc- 
ture^  Topographic  features  cannot  be  seen  since, the  terrain  within 
the  boxed  sub-area  Is  relatively  flat.  However,  the.  displays  show  a 
realistic  depiction  of  a  growing  cumulus  cloud  above  the  corresponding 
topography.  Figure  7-2  shows  initial  development.  Figure  7-6  shows  a 
cumulus  tower  developing  out  of  the  cumulus  clcud. 


Figure  7-1 

Mercator  infrared  image  mS  GMT,  21  June  84 


Figure  7-2 

3-D  cl  Olid- topography  model  for  boxed  sub-area 
in  Figure  7-1. 


Figure  7-3 

Mercator  itifrared  image  for  2229  GMT,  21  June  84. 


Figure  7-4, 

3-D  cloud- topography  model  for  boxed  sub-area 
in  Figure  7-3. 
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Figure  7-5 

Mercator  infrared  image  for  2245  GMT,  21  June  84 


Figure  7-6 

3-D  cloud- topography  model  for  boxed  sub-area 
in  Figure  7-5. 


Figure  7-6  shows  further  development  of  the  cumulus  tower  and  signifi¬ 
cant  development  to  the  east.  The  entire  sequence  shows  the  develop¬ 
ment  of  s  :3mall  isolated  thunderstorm  («<iich  is  propogating  to  the 
east.  Clearing  is  seen  behind  the  storm.  Figures  7-2,  7-4,  and  7-6 
were  part  of  a  four  frame  rovie  loop  which  showed  this  development  in 
3-D  perspective. 

In  cloud  studies,  displays  similar  to  those  shown  in  Figures  7-2, 
7-4,  and  7-6  can  be  used  to  study  vertical  cloud  structure  and  cloud 
development.  Displays  such  as  these  graphically  show  cloud  develop¬ 
ment  in  a  realistic  manner.  Mien  used  in  conjtnction  with  2-D  flat 
satellite  imagery,  greater  image  understanding  is  realized. 

B.'  Briefing  Aids 

Displays  produced  from  the  3-D  cloud-topography  model  contain 
Information  on  cloud  cover  and  weather  systems  which  is  of  interest  to 
pilots.  In  addition,  these  displays  effectively  present  satellite 
data  in  such  a  way  that  the  information  can  be  easily  understood  and 
comprehended.  As  a  result,  these  displays  can  be  used  as  an  aid  to 
brief  pilots  on  weather.  However,  to  make  the  display  more  meaningful 
to  a  pilot  and  to  enhance  the  displays  information  content,  flight 
tracks  can  be  plotted  on  the  3-0  model.  These  flight  tracks  mathema¬ 
tically  define,  in  3-0,  the  path  of  an  aircraft  from  take-off  to 
landing.  Mhen  displayed  with  the  3-0  model,  flight  tracks  enable  the 
pilot  to  visualize  atmospheric  conditions  along  the  planned  route. 

Flight  track;,  can  be  plotted  on  the  3-0  cloud-topography  model 
using  the  FLTRAK  program.  fltRak  accepts  five  parameters  which  are 
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needed  to  mathematically  define  the  3-D  track.  These  parameters 
include  the  aircraft’s  take-off  location,  average  ascent  angle, 
average  cruising  altitude,  average  descent  angle,  and  landing  loca¬ 
tion.  The  take-off  and  landing  locations  are  specified  by  entering 
each  locations  respective  latitude  and  longitude.  Average  ascent 
angle  is  specified  by  entering  the  aircraft’s  clinti  angle  from  the 
point  of  take-off  to  the  point  where  the  aircraft  reaches  its  cruising 
altitude.  Similarly,  average  descent  angle  is  specified  by  entering 
^e  aircraft’s  descent  angle  from  the  point  where  the  aircraft  departs 
from  cruising  altitude  to  the  point  «here  the  aircraft  lands. 

Cruising  altitude  is  specified  by  entering  the  aircraft’s  average 
altitude,  between  ascent  and  descent,  in  meters.  Once  all  parameters 
are  eitered,  rLTRAK  uses  this  data  to  scale  a  model  part  defining  the 
aircraft’s  flight  path.  This  model  part  is  then  merged  to  the  spe¬ 
cified  3-0  cloud-topography  nudel. 

Figure  7-7  shows  a  Mercator  infrared  image  for  2132  GMT,  31 
October  84.  This  figure  shows  the  cloud  cover  for  a  portion  of 
southern  Colorado  and  northern  New  Mexico.  A  line  is  shown  depicting 
a  flight  path  from  Pueblo^  Colorado  (Pl£)  to  Durango,  Colorado  (DRO). 
Figure  7-8  shows  that  this  route  crosses  two  mountain  ranges,  the 
Sangre  de  Cristo  Mountains  to  the  east  and  the  San  Juan  Mountains  to 
the  west.  These  two  ranges  are  separated  by  the  relatively  flat  San 
Luis  Valley.  Figure  7-7  shows  orographically  induced  clouds  formed 
over  these  two  mountain  ranges. 

Figure  7-9  shows  the  corresponding  cloud- topography  display  for 
the  data  shown  in  Figures  7-7  and  7-8.  In  this  figure,  the  3-0 


Figure  7-7 

Mercator  infrared  image  for  3132  GMT,  31  October  84 
with  flight  path  plotted. 


Figure  7-8 

Mercator  topographic  image  wiU  flight 
path  plotted. 
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cloud-topography  model  is  viewed  from  the  southeast.  A  flight  track 
is  shown  depicting,  the  route  from  Pueblo*  Colorado  to  Durango, 

Colorado.  The  track  lias  been  defined  by  a  cruising  altitude  of  6000 
meters  and  ascent  and  descent  angles  of  8  degrees.  In  this  display, 
ascent  and  descent  angles  appear  greater  than  8  degrees  because  the 
vertical  scale  of  the  3-D  cloud-topography  model  has  been  exaggerated 
by  a  factor  of  10.  From  Figure  7-9,  it  is  obvious  that  a' pilot  flying 
this  route  would  encounter  orographi cal ly' induced  clouds  over  the 
mountain  ranges.  However,  at  an  altitude  of  6000  meters,  the  pilot 
should  be  flying  well  above  the  cloud  tops.  In  addition,  the  display 
shows  that  the  pilot  could  expect  clear  conditions  for  take-off  and 
landing.  .  Clear  skies  could  also  be  expected  over  the  San  Luis  valley. 

As  a  briefing  aid,  displays  similar  to  Figure  7-9  could  be  shewn 
to  pilots  to  brief  then  on  weather.  I*"  this  type  of  display, 

topography  data  and  satellite  data  are  presented  in  suc»-  a  way  that 

i 

the  Information  can  be  quickly  and  effectively  comprehended.  Data 
presented  in  this  manner  lets  the  pilot  visualize  corniitions  along  Ihe 
intended  route. 

C.  night  Simulations 

An  alternate  way  of  presenting  cloud  ^d  tojography  data  is 
from  Wie  same  perspective  that  a  pilot  mi^t  see  the  actual  clouds  and 
topography.  This  display  method  uses  KJVIE.BVU  software  to  manipulate 
and  display  the  3-D  model  in  such  a  way  that  simulates  flight- 

Figure  7-10  shows  a  Mercator  visible  image  y*lic  for  2132  GMT,  31 
October  84.  This  image  covers  the  same  geographical  area  shown  in 
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Figures  7-7  and  7-8.  Data  is  also  valid  I'or  f  .e  same  time  period.  A 
line  defining  the  fligiit  path  from  Pueblo.;  Cniorado  to  Durango, 
Colorado  is  also  shown  in  Figure  7-10  h  box  .Jr awn  in  the  center  of 
the  image  defines  a  sub-area  that  is  to  be  trnns formed  int;!  a  3-D 
cloud-topography  model.  Inside  the  boxed-cff  area,  t».o  bands  of 
orographically  induced  clouds  can  be  seen  which  have  formed  over  th-: 
Sangre  de  Cristo  Mountains  to  the  east  and  San  Juan  Mountains  to  the 
west. 

The  3-0  cloud-topography  model  created  for  the  sub-area  in  Figure 
7-10  is  shown  in  Figure  7-11.  Figure  7-11  graphically  shows  the  view 
the  pilot  would  see  ii^iile  flying  at  an  altitude  of  6000  meters  (MSL) 
and  looking  west  along  ^e  line  shown  in  Figure  7-10.  This  view  is 
valid  for  pol*t  A  (labeled  in  Figure  7-10).  In  Figure  7-11,  a  portion 
of  the  Sangre  de  Cristo  Mountains  is  seen  in  the  foreground  and  the 
first. band  of  clouds  formed  above  these  mountains  is  visible.  Off  in 
the  distance,  a  portion  of  the  cloud  free  San  Luis  Valley  is  seen.  On 
the  horizon,  the  secof-d  band  of  clouds  which  are  formed  over  the  San 
Juan  Mountains  is  visible. 

Figure  7-12  shows  the  view  the  pilot  might  see  at  point  B 
.(labeled  in  Figure  7-10).  In  the  foreground,,  cloud  tops  are  visible. 
More  of  the  San  Luis  valley  is  visible  and  the  orographic  clouds  on 
the  horizon  appear  a  bit  closer. 

Figure  7-13  ^ows  the  view  the  pilot  might  see  at  point  C 
(labeled  in  Figure  7-10).  Mere,  the  first  band  of  clwjds  is  no  longer 
visible.  The  cloud  free  San  Luis  Valley  is  seen  in  the  foreground  and 
the  second  band  of  clouds  is  seen  ghi  the  horizon. 


\ 


Figure  7-10 

Mercator  visible  image  for  2132  GMT,  31  October  84 
with  .‘light  path  plotted. 


Figure  7-11 

View  of  3-D  cloud- topography  model  for  boxed 
sub-area  in  Figure  7*10  at  ppint  A. 
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Figure  7-12 

View  of  3-D  cloud- topography  model  for  bo^ed 
sub-area  in  Figure  7-10  at  point  B. 


Figure  7-13 

View  of  3.-0  cloud- topography  model  for  boxed 
sub-area  in  Figure  7-10  at  po'liit  C. 
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The  sequence  (figures  7-11,  7-12,  and  7-13)  graphically  shows 
what  a  pilot  mi^t  see  at  several  points  along  the  hypothetical  flight 
path.  It  must  be  emphasized  that  these  displays  are  created  using 
actual  satellite  data  versus  simulated  data.  Data  displayed  in  this 
way,  or  perhaps  movie  loops  showing  such  a  »quence,  could  be  ^own  to 
pilots  to  inform  them  on  breather  conditions  along  a  fli^t  route. 

The  next  sequence  of  figures  compares  an  actual  photograph  of 
clouds  to  a  computer  graphic  display  of  clouds.  The  photograph  was 


taken  from  a  commerical  airliner  on  approach  to  Stapleton 
International  Airport  in  Denver,  Colorado  on  2  December  84.  The  com¬ 
puter  graphic  display  of  cloud  cover  is  valid  for  2132  GMT,  31  October 
84  for  the  area  shown  in  Figure  6-1.  It  is  inportant  to  note  that  the 
valid  times  of  the  graphic  display  and  the  photograph  do  not  match. 
However,  the  meteorological  situation  ^  similar  in  that  orographically 


induced  clouds  are  formed  over  the  Front  Range  of  the  Rocky  Mountains. 

Figure  7-14  shows  a  scale  model  of  the  cloud  cover  and  topography 
shown  in  Figures  6-1  and  5-1.  The  area  covered  by  this  model  includes 
most  of  central  Colorado  including  the  Denver  area  and  the  Front 
Range. 

Figure  7-15  shows  a  view  of  the  above  cloud-topography  model  from 
the  same  perspective  a  pilot  might  see  clouis  and  topography  on  an 
approa:h  to  Stapleton  International  Airport.  This  view  shows  the 
orograshically  induced  clouds  over  the  Front  Range. 

Figute  7-16  is  a  photograph  of  orographically  irniuced  cloud 
cover  over  ^e  Front  Range.  This  photograph  «»s  takei  looking  west  . 
and  slightly  north,  approximately  100  kilometers  from  the  Front  Range 
using  a  28  millimeter,  wide  angle  lens. 


Figure  7-14 

3-D  cloud-topography  model  for  2132  GMT,  31  October  84 
drawn  to  scale. 
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Figure  7-15  ■ 

3-D  cloud-topography  model  for  2132  GMT,  31  Octobe 
84,  showing  orographically  induced  clouds  over 
the  Frpnt  Range. 


Figure  7-16 

Photograph  of  orographically  induced  clouds  over 
the  Front  Range  taken  from  a  corpnerical  airliner  on 
2  December  84. 


By  comparing  Figures  7-15  and  7-.16,  the  graphic  representation  of 
the  orographic  clouds  is  similar  to  the  actual  photograph  of  the 
orographic  clouds.  This  last  seqience  of  figurei  shows  that  graphi¬ 
cally  produced  cloud  representations  can  be  used  for  fli^t  Simula- 
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Conclusion 

A.  Summary 

In  this  research,  GOES  satellite  data  and  -opography  data 
have  been  mapped  into  high  resolution,  Mercator  images  Mhich  accura¬ 
tely  overlay  each  other.  As  a  result,  cloud  cover  information  ca'.  be 
obtained  with  respect  to  topography,  using  a  matching  set  of  Mercator 
images,  a  sub-area  can  be  specified  which  defines  the  portion  of  the 
images  to  be  transformed  into  a  5-rD,  mathematically  defined  model.  To 
determine  cloud  cover  for  this  sub-area,  a  cloud/no  cloud  threshold  is 
specified  for  the  infrared  image.  Pixels. with  radiance  counts  greater 
than  or  equal  to  the  cloud/no  cloud  threshold  are  considered  cloudy. 
Cloud  top  heights  are  determined  for.  cloudy  pixels  by  n»lating  the 
pixels  radiance  count  to  temperature.  The  resulting  temperatures  are' 
then  related  to  height  using  an  atmospheric  profile  of  temperature 
versus  height.  Cloud  base  heights  are  determined  for  cloudy  pixels  by 
manually  estimating  the  base  height.  This  estimation  is  based  partly 
on  guidance  from  oata  contained  in  the  infrared  image  and  partly  on ' 
meteorological  experience.  Information  contained  in  the  Mercator ! 
topographic  image  is  used  to  produce  a  3-0  topography  model  for  the 
defined  sub-area.  Cloud  top  height  information,  cloud  base  height 


information,  and  cloud  cover  information  is  used  to  produce  3-D  cloud 
nodels  for  the  defined  sub-area.  The  resulting  3-D  topography'  model 
and  the  3-D  cloud  models  are  combined  to  produce  the  3-D  cloud- 
topography  model. 

Three-dimensional  cloud-topography  models  can  be  rotated,  mani¬ 
pulated,  illuminated,  and  displayed  using  MDVIE. BY U  software.  The 
resulting  displays  produce  a  realistic  depiction  of  cloud  cover 
above  a  topographic  surface.  Due  to  the  display’s  visual  appeal  and 
enhanced  information  content,  these  cloud- topography  displays  are  well 
suited  for  cloud  studies,  briefing  aids,  and  flight  simulations. 

For  cloud  studies,  3-0  cloud-topography  models  can  be  used  to 
study  vertical  cloud  structure  and  cloud  development  using  methods 
similar  to  those  currently  used  for  2-D  flat  imagery.  When  used  in 
conjunction  with  2-D  satellite  imagery,  the  displays  produced  from  the 
3-D  model  enable  greater  image  understanding. 

For  pilot  briefing,  3-D  cloud-topography  models  can  be  modified 
to  show  an  aircraft’s  flight  track.  As  a  briefing  aid,  the  displays 
produced  from  the  modified  3-0  model  enable  the  pilot  to  visualize 
atmospheric  conditions  along  the  planned  flight  route.  When  presented 
in  this  manner,  satellite  data  and  topography  data  can  be  quickly  and 
effectively  comprehended. 

For  flight  simulation,  3-0  cluud-topography  models  can  be  rotated 
and  manipulated  in  such  a  way  that  cloud  cover  and  topography  can  te 
viewed  from  the  same  perspective  a  pilot  mi^t  see  the  actual  clcuc 
cover  and  topography.  Data  viewed  in  this  manner  c&n  also  be  used  to 
brief  pilots,  or  could  be  used  in  flight  simulations. 


B.  Future  Work 

0ns  area  for  future  work  would  focus  on  improving  the  3-D 
cloud-topography  model.  Currently,  the  cloud-topogaphy  model  contains 
only  one  layer  of  clouds.  Although  the  software  is  set  up  to  handle 
up  to  three  layers,  current  methods  for  cloud  determination  can  only 
descriminate  one  layer  of  clouds  in  the  satellite  image.  The  addition 
of  the  capability  to,  descriminate  multiple  cloud  layers  would  make  the 
3-D  model  more  realistic.  Another  model  improvement  would  incorporate 
a  quantitative  vertical  scale  into  the  model,  tteckinen  (1984)  descri¬ 
bes  two  excellent  methods  for  scaling  cloud  data  which  could  be 
adapted  to  the  3-b  cloud-topography  model.  The  capabililty  to  display 
a  quantitative  vertical  scale  would  allow  the  user  to  better  estimate 
the  vertical  extent  of  clouds.  Another  model  improvement  would 
incorporate  an  option  to  scale  the  topographic  data  in'  such  a  manner 
that  would  allow  topographic  features  to  be  seen  in  areas  «here 
terrain  is  relatively  flat.  For  example,  the  cloud-topography 
displays  shown  in  Figures  7-2,  7-4,  and  7-6  show  the  development  of  a 
thunderstorm  over  northern  Kansas.  However,  topographic  features  in 
this  area  cannot  be  seen.  The  incorporation  of  advanced  vertical 
scaling  techniques  would  show  topographic  features  without  distorting 
the  clouds  which  are  displayed  above  the  topogr^ic  surface.  One 
final  improve. ent  would  incorporate  more  advanced  techniques  toi  1) 
descriminate  cloudy  and  clear  regions  in  Infrared  satellite  image,  2) 
determine  cloud  top  heig^rts,  and  3)  determine  doud  base  heights. 

These  improvements  would  increase  the  accuracy  of  the  data, being 
displayed. 
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Another  area  for  future  work  would  focus  on  decreasing  the  amount 
of  time  required  to  create  and  displa.'  a  3-D  cloud-topography  model. 
Currently,  it  takes  approximately  five  minutes  to  create  a  3-D  model, 
and  fifteen  minutes  to  display  the  3-D  model.  Thus,  most  of  the  time 
is  spent  displaying  the  3-D  model.  To  significantly  reduce  this  time, 
more  efficient  algorithms  are  needed  to  display  the  3-D  model.  To 
speed  up  the  process,  the  general  purpose  MDVIE.BYU  software  could  be 
abandoned  in  favor  of  now  software  »hich  is  specifically  designed  for 
the  display  of  meteorological  data.  This  new  software  would  use  more 
effici&st  algorithms  to  perform  the  coordinate  transformations  needed 
to  display  the  model.  Another  approach  to  this  problem  would  utilize 
a  nardware/sofware  combination  to  perform  these  coordinate  transfor¬ 
mations.  this  approach  could  resuH  in  systems  that  perform  the  coor¬ 
dinate  transformation  function  and  the  model  display  function  in  real 
time  (approximately  1/30  of  a  second).  A  real  time  system  would  be 
ideal  for  applications  such  as  flight  simulation  »here  actual  movement 
during  flight  could  be  simulated. 

Still  another  interesting  area  for  future  work  would  compare 
actual  photographs  of  clouds  to  comixiter  gr^hic  cloud  representations 
produced  from  the  3-0  cloud-topography  model.  Here,  aircraft  flights 
could  be  made  and  photographs  of  clouds  could  be  taken.  At  tiie  same 
time  the  photographs  are  taken,  satellite  data  could  be  collected. 
Computer  gr^hic  representations  of  clouds  could  be  produced  from 
satellite  data  and  compared  to  the  actual  photogaphs  taken  from  the 
aircraft.  Such  comparisons  would  be  valuable  for  cloud  stiHjies,  model 
improvements,  and  other,  research  artivities. 
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One  last  area  for  future  work  would  explore  additional  applica¬ 
tions  for  the  3-D  cloud-topography  model.  Three  applications  for  the 
3-0  cloud-topography  model  are  discussed  in  this  thesis.  Other  appli¬ 
cations  could  be  developed  *rfiich  would  take  advantage  of  the  MOVIE. BYU 
software  display  capabilities.  One  possible  application  would  incor¬ 
porate  radar  data  into  the  3-0  cloud-topography  model.  Using  this 
idea,  cloud  cover  and  radar  echo  regions  could  be  displayed  to  show 
precipitation  areas  in  clouds.  Another  possible  application  would 
incorporate  other  meteorological  data  into  the  3-D  cloud-topography 
model.  Displays  could  be  produced  showing  cloud  cover  and  isentropic 
surfaces,  cloud  cover  and  pressure  surfaces,  or  cloud  cover  and  any 
other  meteorological  parameter  which  could  be  adequately  described  as 
a  surface.  Another  possible  application  would  use  the  3-D  cloud- 
topography  model  for  spatial  cloud  studies.  Here,  3-D  clouds  could  be 
projected  onto  a  surface  and  percent  cloud  free  line  of  sight  rould  be 
calculated  as  a  function  of  angle  and  distance.  These'  are  but  a  few 
potential  applications  irfiicti  could  be  developed  using  the  3-D  cloud- 
topography  model.  Undoubtedly,  there  are  many  more. 
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